\begin{verbatim}
sage: p = 2^131+1883
sage: fac = (p-1).factor(); fac
2 * 3 * 5 * 37 * 634466267339108669 * 3865430919824322067
sage: primes = [ f[0] for f in (p-1).factor() ]
sage: (p1, p2, p3, p4, p5, p6) = primes
sage: expons = ( (p-1) // r for r in primes )
sage: (n1, n2, n3, n4, n5, n6) = expons
\end{verbatim}%link
